Topological insulators with quaternionic analytic Landau levels 
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We study the 3D topological insulators in the continuum by coupling spin-| fermions to the 
Aharonov-Casher SU(2) gauge field. They exhibit flat Landau levels in which orbital angular 
momentum and spin are coupled with a fixed helicity. Each Landau level contributes one branch 
of gapless helical Dirac modes to the surface spectra, whose topological properties belong to the 
Z2-class. The lowest Landau level states exhibit the quaternionic analyticity as a generalization of 
the complex analyticity of the 2D case. The flat Landau levels can be generalized to an arbitrary 
dimension. Interaction effects and possible experimental realizations are also discussed. 

PACS numbers: 73.43-f,71.70.Ej,73.21.-b 



The 2D quantum Hall (QH) systems p], 0| are among 
the earliest examples of quantum states characterized 
by topology 0, 0] rather than symmetry in condensed 
matter physics. Their magnetic band structures possess 
topological Chcrn numbers defined in time-reversal (TR) 
symmetry breaking systems 0, [Bj]. The consequential 
quantized charge transport originates from chiral edge 
modes [1, 0] , a result from the chirality of Landau level 
wavefunctions. Current studies of TR invariant topo- 
logical insulators (TIs) have made great success in both 
2D and 3D. They are described by a Z2-invariant which 
is topologically stable with respect to TR invariant per- 
turbations [s - 1_7 1 . On open boundaries, they exhibit odd 
numbers of gapless helical edge modes in 2D systems and 
surface Dirac modes in 3D systems. TIs have been experi- 
mentally observed through transport experiments 18- 20j| 
and spectroscopic measurements [2ll427|. 

The current research of 3D TIs has been focusing on 
the Bloch-wave band structures. Nevertheless, LLs pos- 
sess the advantages of the elegant analytic properties and 
flat spectra, both of which have played essential roles in 
the study of 2D integer and fractional QH effects [28j-|36j . 
As pioneered by Zhang and Hu [13], LLs and QH effects 
have been generalized to various high dimensional man- 
ifolds 37H42J. However, to our knowledge, TR invariant 
LLs have not been studied in 3D flat space before. It 
would be interesting to develop the LL counterpart of 
3D TIs in the continuum independent of the band in- 
version mechanism. The analytic properties of 3D LL 
wavefunctions and the flatness of their spectra provide 
an opportunity for further investigation on non-trivial 
interaction effects in 3D topological states. 

In this article, we propose flat LLs in 3D space in 
which spin-i fermions are coupled to an isotropic SU{2) 
Aharonov-Casher potential. When odd number LLs are 
fully filled, the system is a 3D Z 2 TI with TR symmetry. 
Each LL state has the same helicity structure, i.e., the 
relative orientation between orbital angular momentum 
and spin. Just like that the 2D lowest LL (LLL) wave- 
functions in the symmetric gauge are complex analytic 
functions, the 3D LLL ones are mapped into quaternionic 
analytic functions. Different from the 2D case, there is no 



magnetic translational symmetry for the 3D LL Hamil- 
tonian due to the non-Abelian nature of the gauge field. 
Nevertheless, magnetic translations can be applied for 
the Gaussian pocket- like localized eigenstates in the LLL. 
The edge spectra exhibit gapless Dirac modes. Their sta- 
bility against TR invariant perturbations indicates the Z 2 
nature. This scheme can be easily generalize to N dimen- 
sions. Interaction effects and the Laughlin-like wavefunc- 
tions for the 4D case are constructed. Possible realiza- 
tions of the 3D LL system are discussed. 

We begin with the 3D LL Hamiltonian for a spin-i 
non-relativistic particle as 
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where A a a « 



, • ■ ^ab c o"^pr c is a 3D isotropic SU{2) gauge 
with Latin indices run over x, y, z and Greek indices de- 
note spin components f, i; G is a coupling constant and 
ct's are Pauli matrices; V{r) = — ^mwpr 2 is a harmonic 
potential with ujq — \qG\/(2mc) to maintain the flat- 
ness of LLs. A can be viewed as an Aharonov-Casher 
potential associated with a radial electric field linearly 
increasing with r as E(r) x a. H 3D ' LL preserves the TR 
symmetry in contrast to the 2DQH with TR symmetry 
broken. It also gives a 3D non-Abelian generalization of 
the 2D quantum spin Hall Hamiltonian based on Landau 
levels studied in Ref. @. To see the 3D LL structure 
more explicitly, H 3D ' LL can be further expanded as a 
harmonic oscillator with a constant spin-orbit (SO) cou- 
pling as Hzp = w— + ^muiQr 2 =p u>oa ■ L where =F apply to 
the cases of qG > (< 0), respectively. The spectra of 
Eq. [1] in the form of were studied in the context of 
the supersymmetric quantum mechanics (43[. However, 
its connection with Landau levels was not noticed. Eq. 
Q] has also been proposed to describe the electrodynamic 
properties of superconductors 3- 4f|. 

The spectra and eigenstates of Eq. Q] are explained 
as follows. We introduce the helicity number for the 
eigenstate of L ■ a, defined as the sign of its eigen- 
value of the total angular momentum J = L + S, 
which equals ±1 for the sectors of j± — I ± 4, rc- 
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FIG. 1: a) The eigenstates of the 3D harmonic oscillator la- 
beled by total angular momentum j± = I ± i. Following the 
solid diagonal (dashed) lines, these states are reorganized into 
the 3D LL sates with the positive (negative) helicity. b) The 
magnetic translation for the LLL state (I = 0) localized at 
the origin in the case of qG > 0, whose spin is set along an 
arbitrary direction in the a;y-plane. The displacement vec- 
tor 6 lies in the plane perpendicular to spin orientation. The 
resultant state remains in the LLL as a localized Gaussian 
pocket. 

spectively. At qG > 0, the eigenstates are denoted 
as il)n r ;j±,j,;l(r) = Rn r ,i(r)yj ± ,j z -i{^), where the radial 

function is R nr i(r) — r l e "g F(—n r , I + |, ^3-); F is the 

confluent hypergeometric function and Iq = is the 

"magnetic" length; yj±,j z -,i{H) are angular solutions with 
j± = I ± h, respectively. Flat spectra appear with infi- 
nite degeneracy in the sector of j + , where the energy 
dispersion l — (2n r + |)fiwo is independent of I, and 
thus n r serves as the LL index. For the sector of j-, 
the energy disperses with I as E~ l = [2(n r + I) + |]fiwo- 
Similar results apply to the case of qG < 0, where the 
infinite degeneracy occurs in the sector of j_ . These LL 
wavefunctions are the same as those of the 3D harmonic 
oscillator but with different organizations. As illustrated 
in Fig. [T] (a), these eigenstates along each diagonal line 
with the positive (negative) helicity fall into the flat LL 
states for the case of qG > (< 0), respectively. 

Compared to the 2D case, a marked difference is 
that the 3D LL Hamiltonian has no magnetic trans- 
lational symmetry. The non-Abelian field strength 
grows quadratically with r as Fij(f) = diAj — djAi — 
H [Ai , Aj] = ge ljk {a k + -^-r k (a-r)}. Nevertheless, mag- 
netic translations still apply to the highest weight states 
of the total angular momentum J = L + S in the LLL at 
qG > 0. For simplicity, we drop the normalization factors 
of wavefunctions below. For the positive helicity states 
with j z = L and 5* are parallel to each other. Their 

wavefunctions are denoted by ip^i(r) — (x + iy) l e u g <gi 

where cXq is the spin eigenstate of ft ■ a with eigen- 
value 1. For these states, the magnetic translation is 
defined as usual T^(5) = exp[— S ■ V + jp-r xy ■ (z x 5)], 

where S is the displacement vector in the xy-plane and 



r xy is the projection of r in the xy-plane. The resultant 

state, T s (5)ip%y(f) = e u l ip^(f-S), remains in the 
LLL. Generally speaking, the highest weight states can 
be defined in a plane spanned by two orthogonal unit 

vectors as i^^{r) = [(ei + ie 2 ) ■ f] l e "g ® ag 3 with 
63 = ei x 62. The magnetic translation for such states 
is defined as Tg 3 (S) = exp[— S ■ V + jp- f\2 ■ {e% x $)], 

where S lies in the ei,2-plane and f\2 — f — e%{r ■ £3). 
However, for other LLL states with j z ^ ±j at qG > 
and all the LLL states at qG < 0, their spin and orbit 
angular momenta are entangled. The above definition of 
magnetic translation fails. As an example, let us trans- 
late the LLL state localized at the origin as illustrated 
in Fig. [T] (b). We set the spin direction of Vj^&o m 
the zy-plane parameterized by £3(7) = x cos 7 + ysin^f, 
i.e., Qe 3 (7) = t) + e* 7 | |)), and translate it along 
e\ = z at the distance R. The resultant states read 

as 1p 7tR (p,<p,z) = e 'fflpsin(0- 7 ) e -|r-fl£| 2 /44 (g)Q!g 3 ( 7 ), 

where p = y/x 2 + y 2 and <j> is the azimuthal angular of r 
in the xy-plane. Such a state remains in the LLL as an 
off-centered Gaussian wave packet. 

The highest weight states and their descendent states 
from magnetic translations defined above have clear clas- 
sic pictures. The classic equations of motion are derived 
as 

r = — p + 2oj (rx —S), p — 2u>opx —S — muiZr, 
m n n 

S=^SxL, (2) 

where p is the canonical momentum, L = r x pis the 
canonical orbital angular momentum, and S here is the 
expectation value of ^<r. The first two describe the mo- 
tion in a non-inertial frame subject to the angular veloc- 
ity ■^■S, and the third equation is the Larmor precession. 

L ■ S is a constant of motion of Eq. [2] In the case of S \\ L, 
we can prove that both S and L are conserved. Then the 
cyclotron motions become coplanar within the equatorial 
plane perpendicular to S. Centers of the circular orbitals 
can be located at any points in the plane. 

The above off-centered LLL states break all the ro- 
tational symmetries. Nevertheless, we can recover the 
rotational symmetry around the axis from the origin to 
the packet center by performing the Fourier transform 
over ip lt n(p, cj), z) with respect to the azimuthal angle 7 
of spin polarization. The resultant state is a j 2 -eigenstate 

cx e^f^e™*{ J m (x)\ t) + J„ i+ i(x)e^| |)}, (3) 

with x = Rp/(21q). In fact, ^ 7j r's are non-orthogonal 
to each other, thus this transformation is linear but not 
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unitary. For R > 0, although Eq. [3] exists for all the 
values of m, the norm of the direct result of this trans- 
formation decays as increasing \m\. At large distance of 
R, the spatial extension of ipj z=m+ i ji in the xy-plane 
is at the order of uiIq/R, which is suppressed at large 
values of R and scales linear with m. In particular, the 
narrowest states ip±i- r exhibit an ellipsoid shape with 
an aspect ratio decaying as Ig/R when R goes large. For 
those states with \m\ < R/Iq, they localize within the 
distance of Iq from the center Rz. As a result, the real 
space local density of states of LLL grow linearly with R. 
In fact, this local density of state can be calculated ex- 



as r —> - 



(- 



Ig 



2r i 



which 



actly as p(r) 

approaches as r — > +oo. 

In analogy to the fact that the 2D LLL states are com- 
plex analytic functions due to chirality, we have found 
an impressive result that the helicity in 3D LL systems 
leads to the quaternionic analyticity. Quaternion is the 
first discovered non-commutative division algebra, which 
has three anti-commuting imaginary units i, j and k, sat- 
isfying i 2 = j 2 = k 2 = — 1 and ij = k. It has been applied 
in quantum systems 0, EH and SO coupled BECs (49j . 
Just like two real numbers forming a complex number, 
a two-component complex spinor ip = (ipf,ifii) T can be 
viewed as a quaternion defined as f — ifo + jtpp In the 
quaternion representation, the TR transformation ia2ip* 
becomes Tf = —fj satisfying T 2 = — 1; multiplying a 
U(l) phase factor e^ip corresponds to /e 1 *; the SU(2) 
operations e~ l ~5~^i\), e _l ^ i ^'0j an d e -1 ^"^^ map to e?*/, 
e^/, and e~i^f, respectively. The quaternion version 
LLL (x,y,z) = >!'•., ., • .j • 

1 r'Yi,. . 

which is a 

* ■ - / / 

quaternionic polynomial. 

As a generalization of the Cauchy-Riemann condition, 
a quaternionic analytic function f(x,y,z,u) satisfies the 
Fueter condition [5Cj| as 



oiib LLL . 



where 



(j+,j z \l,m 



1:1 

' 2' 2- -7— ■ - 

Gaussian factor does not appear in /f^f 



2 J 21' x l,j+,3* ~ 

Please note that the 
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(4) 



where x, y, z and u are coordinates in the 4D space. In 
Eq. |U imaginary units are multiplied from the left, thus 
it is the left- analyticity condition which works in our con- 
vention. Below, we prove the LLL function (x, y, z) 
satisfying Eq. 01 Since is defined in 3D space, it 

is a constant over u, and thus only the first three terms 
in Eq. 2] apply to it. Obviously the highest weight states 
with spin along the z-axis, fj—j — j+i = ( x + iy) > sat- 
isfy Eq. U which is reduced to complex analyticity. By 
applying an arbitrary SU(2) rotation g characterized by 
the Eulerian angles (a, /3, 7), fj^=j z transforms to 



1, LLL 



(x,y,z) 



= p. '2 



e^e-^f^(x',y\z% (5) 



where (x',y', z') are the coordinates by applying the 
inverse of g on (x,y,z). We check that + 



e 2 e J 2 e 2 



/_8_ 

1 dx' ^ 6 dy' ^ 

i^}fj+,fS x ' , >y' , > z ') ~ ®' Essentially, we have proved 
that a quaternionic analytic function remains analytic 
under rotations. Since all the highest weight states are 
connected through SU(2) rotations, and they form over- 
complete basis for the angular momentum representa- 
tions, we conclude that all the 3D LLL states with the 
positive helicity are quaternionic analytic. 

Next we prove that the set of quaternionic LLL states 

ff+=i+i-,jSi+ - 5' Jz - J '+ - form thc com P lete 

basis for quaternionic valued analytic polynomials in 3D. 
We only need to prove it for each value of j+ = I + \. 
Any linear superposition of the LLL states with j + can 

be represented as ft = Ysjt=-j f-f+.fz c j* > wnere c j z is a 
complex coefficient. Because of the TR relation /; 

pLLL 
j+i-. 

independent basis as 



LLL _ 

ff+,-j m h fi can be expressed in terms of / + 1 linearly 



fi(x,y,z) 



EfLLL 



(6) 



where q m — c m+ 1 —jc_ m _i is a quaternion constant. On 
the other hand, it can be calculated that the rank of the 
linearly independent l-th order quaternionic polynomi- 



als satisfying Eq. [4] is just 



J+l, 



thus fj^jl 's with j z > i are complete. 

The topological nature of the 3D LL problem exhibits 
clearly in the gapless edge states. Inside the bulk, LL 
spectra are flat with respect to I. As I goes large, the 
classical orbital radius r c approaches the boundary with 
the radius i?o- For example, for a LLL state, r c = \/~211g. 
As I > l c 1=3 ^(Ro/Ig) 2 , the wavefunctions concentrate on 
the boundary and become surface states. Their spectra 



become E(l) « 1(1 + 1) 



•ImR't 



— IhujQ. When the chemi- 



cal potential /i lies inside the gap, it meets the surface 
states with the Fermi angular momentum denoted by If. 
The edge spectra can be linearized around If as Hm — 
(v f /R )a -L-fi with v f = (21 f + l)h/(2MR ) - ojR . 
This is the Dirac equation defined on a sphere with the 
radius Rq. It can be expanded around r = Roi r as 
Hbd = tiVf(k x (?) • e r — /i. Similar reasoning applies 
to other Landau levels which also give rise to Dirac spec- 
tra. Due to the lack of Bloch wave band structure, it 
remains a challenging problem to directly calculate the 
bulk topological index. Nevertheless, the Z 2 structure 
manifests through the surface Dirac spectra. Since each 
fully occupied Landau level contributes one Dirac cone 
on the surface, the bulk is Z2-nontrivial if odd number of 
LLs are occupied. 

The above scheme can be easily generalized to arbi- 
trary dimensions by combining the 7V-D harmonic oscil- 



4 



lator potential and SO coupling. For example, in 4D, we 



r 4D 



W oEl<a<ii<4 r ° i ' i « 



have H iD > LL = ^ + f mwgr; 
where L a b = r a pb~rbp a and the 4D spin operators are de- 
fined as r« = -\[a\a% T li = ±a l with 1 < i < j < 3. 
The ± signs of T z4 correspond to two complex conjugate 
irreducible fundamental spinor representations of 50(4), 
and the + sign will be taken below. The spectra of the 
positive helicity states are flat as E +yTlr — (2n r + 2)huj. 
Following a similar method in 3D, we prove that the 
quaternionic version of the 4D LLL wavefunctions sat- 
isfy the full equation of Eq. 2) They form the com- 
plete basis for quaternionic left-analytic polynomials in 
4D. For each l-th order, the rank can be calculated as 

('t 3 )-(T)^G + i)G + 2). 

We consider the interaction effects in the LLLs. For 
simplicity, let us consider the 4D system and the short- 
range interactions. Fermions can develop spontaneous 
spin polarization to minimize the interaction energy in 
the LLL flat band. Without loss of generality, we as- 
sume that spin takes the eigenstate of T 12 = T 34 = <r 3 
with the eigenvalue 1. The LLL wavefunctions satisfy- 
ing this spin polarization can be expressed as = 

(x + iy) m {z + iu) n e~^ (gi | a ) with |a) = (1,0) T . The 4D 
orbital angular momentum number for the orbital wave- 
function is / = m + n with to > and n > 0. It is easy to 
check that > 4D is the eigenstate of £ afc L ab T ab with 

the eigenvalue (m + n)K. If all the :4D, s are filled 
with < to < N m and < n < N n , we write down a 
Slater-determinant wavefunction as 



*(«i,ioi; • • • ; v N , w N ) = det[vi 



(7) 



where the coordinates of the i-th particle form two pairs 
of complex numbers abbreviated as i>j = xi + iyi and 
Wi = Zi+iui] a, P and i satisfy < a < N m , < /3 < N n , 
and 1 < i < N — N m N n . Such a state has a 4D uniform 
density as p — -rrn ■ We can write down a Laughlin-likc 
wavefunction as the fc-th power of Eq. [7]whose filling rel- 
ative to p should be 1/k 2 . For the 3D case, we also con- 
sider the spin polarized interacting wavefunctions. How- 
ever, it corresponds to that fermions concentrate to the 
highest weight states in the equatorial plane perpendic- 
ular to the spin polarization, and thus reduces to the 
2D Laughlin states. In both 3D and 4D cases, fermion 
spin polarizations are spontaneous, thus low energy spin 
waves should appear as low energy excitations. Due to 
the SO coupled nature, spin fluctuations couple to orbital 
motions, which leads to SO coupled excitations and will 
be studied in a later publication. 

One possible experimental realization for the 3D LL 
system is the strained semiconductors. The strain ten- 
sor e a b = \{d a Ub + dbUa) generates SO coupling as 

HsO = fo&[(£xyky~£xzkz)&x + {€zyk z — €xykx)c r y-\~{£zxkx~ 

£yzky)<Jz\ where a = 8 x 10 5 m/s for GaAs. The 3D strain 



configuration with u = ^(yz, zx,xy) combined with a 
suitable scalar potential gives rise to Eq. [T]with the cor- 
respondence wo = \ a i '■ A similar method was proposed 
in Ref. Q to realize 2D quantum spin Hall LLs. A LL 
gap of lmK corresponds to a strain gradient of 1% over 60 
/im, which is accessible in experiments. Another possible 
system is the ultra-cold atom system. For example, re- 
cently evidence of fractionally filled 2D LLs with bosons 
has been reported in rotating systems |51| . Furthermore, 
synthetic SO coupling generated through atom-light in- 
teractions has become a major research direction in ultra- 
cold atom system [HD, [53[. The SO coupling term in 
the 3D LL Hamiltonian uja ■ L is equivalent to the spin- 
dependent Coriolis forces from spin-dependent rotations, 
i.e., different spin eigenstates along ±x, ±y and ±z axes 
feel angular velocities parallel to these axes, respectively. 
An experimental proposal to realize such an SO coupling 
has been designed and will be reported in a later publi- 
cation [HH- 

In conclusion, we have generalized the flat LLs to 3D 
and 4D flat spaces, which are high dimensional topo- 
logical insulators in the continuum without Bloch-wave 
band structures. The 3D and 4D LLL wavefunctions 
in the quaternionic version form the complete bases of 
the quaternionic analytic polynomials. Each filled LL 
contributes one helical Dirac Fermi surface on the open 
boundary. The spin polarized Laughlin-like wavefunc- 
tion is constructed for the 4D case. Interaction effects 
and topological excitations inside the LLLs in high di- 
mensions would be interesting for further investigation. 
In particular, we expect that the quaternionic analyticity 
would greatly facilitate this study. 

This work grew out of collaborations with J. E. Hirsch, 
to whom we are especially grateful. Y. L. and C. W. 
thank Y. S. Wu for helpful information on quaternions, 
S. C. Zhang and J. P. Hu for helpful discussions. Y. 
L. and C. W. are supported by NSF DMR-1105945 and 
AFOSR. FA9550-11-0067 (YIP). 

Note Added: Near the completion of this manuscript, 
we learned that the 3D Landau level problem is also stud- 
ied by S. C. Zhang (5^. 
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FIG. 2: The algebra structure of the 3D Landau levels in 
the positive helicity sector. Operators A± connect states with 
different I in the same Landau level, while B- and C+ connect 
those between neighboring Landau levels. 
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Supplementary material for 
"Topological insulators with 
quaternionic analytic Landau levels" 

In this supplementary material, we present several 
points including the ladder algebra to explain the de- 
generacy of the 3D Landau level (LL) wavefunctions, the 
quaternionic version of the 3D lowest Landau level (LLL) 
states with the negative helicity, the numerical calcula- 
tion on the 3D LL spectra with the open boundary, and 
the generalization of LLs to an arbitrary dimension. 

We have shown in the main text that there is no mag- 
netic translation symmetry for 3D LL Hamiltonian. Nev- 
ertheless, the spectra flatness can still be understood by 
constructing the following ladder algebra structures. For 
example, we take the case of qG > and consider the 
positive helicity Landau level states of H + . The vari- 
able transformation for the radial eigenstates is applied as 
Xn,.,;( r ) = r -Rn,.,;( r )j and the corresponding radial Hamil- 
tonians become 



Hi = hw (- 



1(1 + 1) 



dr* 



(8) 



where the dimensionless radius reads as r* = t/Iq. The 
ladder operators are defined as 



1 * 

r 



. d l + l 

v ' dr* r* 2 



(9) 



They satisfy the relations -H"z±i^4±(0 = A±(l)Hi. Con- 
sequentially, Xn r ,i±i = A±(l)xn r ,l w hh the same energy 
independent of I. All the states in the same LL can be 
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FIG. 3: The energy dispersion of the first four Landau levels 
v.s. I = j — Open boundary condition is used for a ball with 
the radius Ro/Ig = 8. The edge states correspond to those 
with large values of I and develop linear dispersions with I. 
The most probable radius of the LLL state with I is r — IgVI- 

reached by successively applying A± operators. To con- 
nect different LLs, other two ladder operators are defined 
as 



c +<" = ^-^ + r*' < 10 » 

which satisfy = B_(l)(Hi + 2huj ) and 

Hi + iC+(l) = C+(l)(Hi — 2frw ), respectively. By apply- 
ing B-(l) (C+(l)) to X„,.,/(r), we arrive at Xn r +i,i-i = 
B-(l)Xn r ,i and Xn r -i,i+i = C + (l)x nr ,l witn the energy 
change ±2fi,wo, respectively, as illustrated in Fig. [5] Simi- 
lar algebra can also be constructed for the case of qG < 0. 

We have numerically calculated the spectra of the 3D 
LL Hamiltonian with the open boundary condition for 
the positive helicity states with j + = £ + i. The results 
for the first four LLs are plotted in Fig. [5] The radius 
of the boundary is Ro/Ig = 8. For the lowest LL (LLL) 
states, when the orbital angular momentum I exceeds a 
characteristic value l c w 30, the spectra become disper- 
sive indicating the onset of surface states. 

In the main text, we have showed that the 3D LLL 
states with the positive helicity in the quaternion repre- 
sentation form a set of complete basis for the quaternionic 
left-analytic polynomials. For the case of the LLL with 
negative helicity, their quaternionic version gf_f {%> y, z) 
are not analytic any more. Nevertheless, they are related 
to the analytic one through Q^/tn^si ^ — (~xk — yj + 

The study in 3D and 4D LL systems can be gener- 
alized to N-D by replacing the vector and scalar po- 
tentials in Eq. 1 in the main text with the SO(N) 
gauge field A a (f) = gr b S ab and V{r) = -^mLu r 2 , 
respectively, where S ab are the SO(N) spin operators 
constructed based on the Clifford algebra. The rank-fc 
Clifford algebra contains 2k + 1 matrices with the di- 
mension 2 k x 2 k which anti-commute with each other 
denoted as T" (1 < a < 2k + 1). Their commutators 
generate r ab = -^[T a ,T b ] for 1 < a < b < 2k + 1. 



For odd dimensions N = 2k + 1, the SO(N) spin op- 
erators in the fundamental spinor representation can be 
constructed by using the rank-fc matrices as S ab = ^T ab . 
For even dimensions N = 2k + 2, we can select 2fc + 2 
ones among the 2fc + 3 T-matriccs of rank-(fc + 1) to 
form S ab = ir ab , then all of S ab commute with T 2k + 3 . 
This 2 fc+1 -D spinor Rep. of S ab is thus reducible into 
the fundamental and anti-fundamental Reps. Both of 
them are 2 fc -D, which can be constructed from the rank- 
fc r-matrices as S°" 2k+2 = ±\T a (l < a < 2k + 1) and 
S ab _ l r a b (l < a < b < 2k + 1), respectively. As 
for TR properties, r a 's are TR even and odd at even 
and odd values of fc, respectively. We conclude that at 
N = 2k + 1, the 7V-D version of the LL Hamiltonian is 
TR invariant in the fundamental spinor Rep. At N = 4fc, 
it is also TR invariant in both the fundamental and anti- 
fundamental representations. However N = 4fc + 2, each 
one of the fundamental and anti-fundamental Reps is not 
TR invariant, but transforms into each other under TR 
operation. 

Similarly, the A^-D LL Hamiltonian can be reorganized 
as the harmonic oscillator with SO coupling. For the case 
of qG > 0, it becomes 

H N ,+ = + ^muj 2 r 2 - Hoj T ab L ab , (11) 

where L a b = r a pb — rtp a with 1 < a < b < N. The l-th 
order A^-D spherical harmonic functions are eigenstates 
of L 2 = L a i,L a b with the eigenvalue of h 2 l(l + D — 2). 
The N-D harmonic oscillator has the energy spectra of 
E nrt i = (2n r + 1 + N/2)hu. When coupling to the funda- 
mental spinors, the l-th spherical harmonics split into the 
positive helicity and negative helicity (j_) sectors, 
whose eigenvalues of the T ab L ab are HI and — h(l + N — 2), 
respectively. For the positive helicity sector, its spectra 
become independent of I as E+ = (2n r + Nj2)fko, with 
the radial wave functions are 

RnAr) = r l e- r2 / 4l oF(~n r , I + N/2, r 2 /2l 2 G ). (12) 

The highest weight states in the LLL can be written as 

iPabM^ = ± ■ rf e-*- 2 / 4 ^ ® a ± , lbl (13) 

where a±, ab is the eigenstate of T a b with eigenvalue ±1, 
respectively. The magnetic translation in the a6-plane 
by the displacement vector S takes the form T ab (S) = 
cxp[— 5 ■ V + 2p~T a b(r a Sb — TbSa)]- Similarly to the 3D 
case, staring from the LLL state localized around the ori- 
gin with I = 0, we can perform the magnetic translation 
and Fourier transformation with respect to the transverse 
spin polarization. The resultant localized Gaussian pock- 
ets are LLL states of the eigenstates of the SO(N — 1) 
symmetry with respect to the translation direction i5. 
Again each LL contributes to one channel of surface Dirac 
modes on S ,JV ~ 1 described by Hbd = (vf /Ro)T a bL a b — [i. 



